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EQUATORIAL RADIUS OF THE EARTH: 
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ABSTRACT 

An interesting variation on the familiar method of determining 
the earth's equatorial radius a e , from a knowledge of the earth's 
equatorial gravity is suggested. The value of equatorial radius 
thus found is 6378,142 ± 5 meters. The associated parameters 
are GM = 3.986008 ± .000004 x 10 20 cm 3 sec" 2 which includes 
the relative mass of atmosphere yl0~ 6 xGM, the equatorial gravity 
y e = 978,030.9 milligals (constrained in this solution by the 
Potsdam Correction of 13.67 milligals as the Potsdam Correc- 
tion is more directly, or less indirectly , measurable than the 
equatorial gravity) and an ellipsoidal flattening of f = 1/298.255. 
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EQUATORIAL RADIUS OF THE EARTH: 

A DYNAMICAL DETERMINATION .. 

INTRODUCTION 

The problem of determining the earth's equatorial radius _a , consists of scaling 
a mean ellipsoid which best fits an accurately determined detailed geoid. The 
importance of the problem is well known in geodesy and many of the geodetic 
and geophysical applications of artificial earth satellites. Several theoretical 
approaches are available to solve the problem, namely, gravimetric methods 
(see for example, Heiskenan and Moritz, 1967), astrogeodetic methods (see for 
example, Fricke et al., 1965), dynamic and geometric satellite methods (see for 
example, Kaula, 1967; Anderle, 1967 or Veis, 1967) but their applications is 
limited by the density and precision of the available data required for the purpose. 

In presatellite era a value of a e = 6378,388 meters, computed by Hayford (1909) 
from the isostatically reduced astrogeodetic data in the United States, was 
adopted to define the equatorial radius of the "international. reference ellipsoid." 
The artificial earth satellite data analyses and space probe studies have indicated 
much smaller values (see, for example, Veis, 1967) ranging approximately be- 
tween a = 6378,124 meters and a = 6378,160 meters, with the latest determi- 
nations concentrated around a value of a g = 6378,140 meters approximately. 

This can be seen from Table 1 which, by the way, is by no means anywhere near 
a complete listing of all the important work on the determination of a e . 
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Table 1 

Comparisons of the Recent Estimates 
of the Earth's Equatorial Radius a^ 


Reference 

Equatorial Radius 
a e (meters) 

Vies (1965) 

6378,169 

Geodetic Reference System, 1967 

6378,160 

Kaula, 1967 

6378,153 

Present paper (atmospheric mass included) 

6378,144 

Present paper (atmospheric mass excluded) 

6378,142 

Veis (1965) 

6378,142 

Lambeck (1971) 

6378,140 

Rapp (1970) 

6378,138 

Strange et al. (1971) 

6378,124 


THE METHOD 

The method suggested here stems from Khan's (1972) work showing that the 
Potsdam Correction d, can be expressed as 

d = (Tp 1 - 7 s ) + (Ag 1 - Ag s ) (1) 

The quantities y and Ag are the theoretical gravity and gravity anomaly on the 
appropriate reference surface respectively. The superscript I denotes the 
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international reference ellipsoid and the superscript S the appropriate "satellite 
reference ellipsoid." The subscript P refers to the "Potsdam' Gravity System." 
The theoretical gravity y on the surface of an arbitrary reference ellipsoid, is 
given by the well-known formula 

7 ~ 7 e ( 1 + a 2 sin 2 4> + a 4 sin 2 2 4>) (2) 

where the equatorial gravity y e is related to GM, a e and the flattening f as 


GM = y e a 2 (1 “ f) 1 + - m + ■= m f + 


9. 


m 


and the coefficients a 2 and a 4 are given as 


5 1 2 26 r 15 , 

a- 2 —f + j m + 2 ' 7 mf + J m 2 


and 


1 r9 5 

a 4 = g f ‘ g m f 


with 


w 


m 


2 a 3 (1 ~ f) 
GM 


(3) 


(4) 


Note that the above formulas are accurate only to the O (f 2 ) as they neglect , 
quantities 0(f 3 ). 

The appropriate choice of y e , a 2 and a 4 in the above equations will yield y 1 or 
y s as the. case may be. In practice, the quantity 7p is computed from the 
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international gravity formula in which the value of equatorial gravity is tied 
to the Potsdam gravity system. Hence any error in the system is reflected 
in this quantity. The quantity y s is computed from the socalled "satellite gravity 
formula" (not the "Gravity Formula 1967") in which the ellipsoid flattening is 
held fixed at f = 1/298.255 but GM and a e are treated as variables. 

The gravity anomaly Ag is given by 


a * n= 2 m=0 


n + 2 


GM 00 n / a e 

A g 7 £ £ It) ( n _ !) ( S C nra “S mk + S S rm sin m M L( sin <£) (5) 


where 
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g s - observed g - reference g 


and other symbols have the usual meaning (see for example, Khan, 1969). In 
practice, however, only C 20 and C 40 are considered for the reference ellipsoid. 
Consequently all other S C nm and § S nm (except n = 2 and 4) equal the observed 
geopotential coefficients. The largest of the coefficient differences, involved in 
the computation of A g 1 and A g s or the difference Ag 1 - Ag s , does not exceed 
10 _s . Hence, neither of these quantities is sensitive to small changes in GM 
and a e . Thus the sensitivity of 7 s to small variations in GM or a e or more 
exactly the ratio GM/a e 2 , on which the equatorial gravity y e s , used in computing 
7 s , is based makes the relation stated in Equation (1) highly useful as it can be 
utilized eventually to determine a e . 

Note from Equation (3) and the expression for m that the equatorial gravity y e is 
automatically fixed once the earth's rotational speed co, GM, a e and the flattening 
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f are defined. Of these quantities w and f are well-determined. Of the remaining 
three i.e., GM, a e and y e , knowledge of any two will determine the third uniquely. 
The quantity GM is independently determined from lunar and planetary probes 
and although these determinations still have a considerable uncertainty asso- 
ciated with them, GM is far more strongly determined than y e or a e . The 
quantity y e is not directly measurable and must be either derived theoretically 
from expressions involving a e (see Equation 3) which we wish to compute in this 
paper or from surface gravity data which does not form an adequate network. 

(This was true until recently but the situation is likely to change with the increased 
use of the portable free fall absolute gravity measuring. equipment developed by 
Hammond and Faller, 1971). 

This difficulty can be overcome, however, by controlling the value of equatorial 
gravity y e s (entering Equation (1) via 7 s ) through the Potsdam gravity correction 
d as demonstrated in Equation (1) . Because of the recent advances in electronic 

technology, the absolute gravity value at a particular site can be measured to 
an accuracy of a few parts in 10 8 of absolute gravity g (Hammond and Faller, 
1971). If such precise measurements could be made at the Potsdam site itself, 
the Potsdam gravity correction could be determined with the same precision. 
However, since no such measurements are available (in case they are, they have 
not been released to the scientific community), we can use indirect estimates of 
d made through highly accurate measurements of absolute gravity at sites pre- 
viously-tied to the Potsdam gravity system. Such measurements yield (Cook 
and Hammond, 1969) a value of approximately d = 13.67 milligals which is 
probably correct to better than a milligal. 
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The advantage of the method suggested in Equation (1) lies in the fact that the 
less cooperative equatorial gravity value is controlled through the Potsdam 
Correction d which is relatively more directly, or in fact less indir ectly, 
measured. This is clearly shown by Equations (1) through(4) . These equations 
also demonstrate that the method is equivalent to applying the Potsdam Correc- 
tion to the equatorial gravity value y* in the international gravity formula, 
projecting the corrected value onto an ellipsoid with a flattening of f = 1/298.255 
and using the resulting value as a known parameter in Equation (3) to determine 
a e . Thus Equation (1) is really an interesting variation on the well-known rela- 
tion stated in Equation (3). 

RESULTS 

The solution of Equation (1) is given in Figure 1 which shows plots of the earth's 
equatorial radius a e as functions of GM and Potsdam gravity correction d. The 
quantities y p x and Ag 1 are based on the well-known parameters of the interna- 
tional gravity formula and the international reference ellipsoid. The quantities 
y s and Ag s are based on an ellipsoidal flattening of f s = 1/298.255. The value 
of Ag 1 - Ag s is kept compatible with the indicated values of GM and a e , though 
small discrepancies in the value of GM and a e will not affect Ag 1 - Ag s to the 
degree of accuracy (a few parts in 10 8 of absolute gravity g) considered here. 

All the parameters used in computing y s , which makes chief contribution to the 
changes in d (see Equation (1) and Figure 1), are computed exactly from the 
appropriate values of GM and a e . 

As pointed out earlier, a well-accepted value of Potsdam correction (Cook and 
Hammond, 1969; Khan, 1972), obtained from the free fall determinations of 
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GRAVITATIONAL CONSTANT * MASS OF THE EARTH, GM 
no 20 CM 3 SEC- 2 ) 






absolute gravity acceleration at various sites, is 13.67 milligals, with an accuracy 
probably better than a milligal. The most recent value of GM is 3.986008 x 
10 20 cm 3 sec -2 (Esposito and Wong, 1972). Corresponding to these quantities, 
the value of equatorial radius is found to be a e = 6378,144 meters (to the nearest 
meter). The correction for the mass of the atmosphere reduces GM to 3.986005 
x 10 20 cm 3 sec -2 and the corresponding value of the equatorial radius is a e = 
6378,142 meters (to the nearest meter). Compare it with a e = 6378,160 meters 
adopted by IUGG (Moscow, August 1971) on one extreme and with a e = 6378,124 
meters (Strange et al., 1971) on the other extreme. Some of the important recent 
values of a e are listed in Table 1. 

The maximum error introduced by the approximations made in the second order 
theory used here is 5 x 10 -8 milligals. The error contributed by this factor in 
the determination of a e , therefore, cannot exceed a few tenths of a meter. The 
chief sources of error are the uncertainties in GM, Potsdam Correction and the 
value of equatorial gravity in the International Reference System. Of these the 
uncertainty in GM is ±4 x 10 14 cm 3 sec -2 . This corresponds to an uncertainty 
of ±3.2 meters in the value of a e . The maximum uncertainty in Potsdam Cor- 
rection is not believed to exceed 1 milligal. This corresponds to a maximum 
uncertainty of 3 meters. Additional uncertainties could stem from some yet 
undetected systematic shift in the Potsdam gravity correction or the value of 
absolute gravity acceleration at the Potsdam site but it seems unlikely as the 
Potsdam correction obtained from various observation sites remains consistent. 
The tacit assumption made here is that the gravity value obtained from the in- 
ternational gravity formula at the projection of the Potsdam site onto the 
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reference ellipsoid is compatible with the old gravity value at the Potsdam site 
i.e., no error is assumed in the value of equatorial gravity yj- derived by 
Heiskenan (1928), except for the Potsdam Correction, and in the reduction of 
the absolute gravity value at the Potsdam site to the geoid. With this granted, 
the value of the equatorial radius of the best-fit earth ellipsoid is found to be 
a = 6378,144 ± 5 meters with the mass of the atmosphere included in GM and 

e 

a e = 6378,142 ± 5 meters with GM corrected for the mass of the atmosphere. 
The corresponding value of equatorial gravity y e compatible with a e = 6378,142 
meters is y e = 978,030.9 milligals (rounded to the nearest tenth of a milligal). 
Compare it to the equatorial gravity value of 978,031.1 milligals (see for ex- 
ample, Veis, 1967) obtained from a much more extensive gravity coverage than 
was available to Heiskenan (1928). It is remarkable that y e should display such 
a stability. It is also remarkable that if Heiskenan had the present values of 
GM, f and d available to him in 1928, he would have obtained a value of a e very 
close to the modern values of the earth's equatorial radius. 

RECOMMENDATIONS 

The Recommended Geodetic Parameters are: 


GM (corrected for mass 
of the atmosphere) 


3.986005 ± .000004 x 10 2Q cm 3 sec -2 
6378,142 ± 5 meters 
298.255 ± .005 


d = 13.67 ± 1 milligals 
y e = 978,030.9 ± 2.5 milligals 
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LIST OF ILLUSTRATIONS 


Figure 1. Earth's equatorial radius a e as function of GM and d. The labels on 
a e curves should be read with 6378 as prefix. For example, the 
rightmost slant curve has a value of 6378,161 meters and the leftmost 
6378,136 meters. 
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